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Consider a random medium consisting of points randomly distributed so that there is no corre- 
lation among the distances. This is the random link model, which is the high dimensionality limit 
(mean field approximation) for the euclidean random point structure. In the random link model, 
at discrete time steps, the walker moves to the nearest site, which has not been visited in the last 
H steps (memory), producing a deterministic partially self avoiding walk (the tourist walk). We 
have obtained analitically the distribution of the number n of points explored by a walker with 
memory = 2, as well as the transient and period joint distribution. This result enables to explain 
the abrupt change in the exploratory behavior between the cases j-i — 1 (memoryless, driven by 
extremal statistics) and ^ = 2 (with memory, driven by combinatorial statistics). In the /i = 1 case, 
the mean newly visited points in the thermodynamic limit (A*' ^ 1) is just (n) = e = 2.72 . . . while 
in the (1 = 2 case, the mean number (n) of visited points is proportional to N^^^ . Also, this result 
allows us to stabilish an equivalence between the random link model with /i = 2 and random map 
(uncorrelated back and forth distances) with /i = and the drastic change between the cases where 
the transient time is null compared to non-null transient times. 



I. INTRODUCTION 

Although not as thoroughly studied as random walks 
in disordered media [l| and complex media 0, which 
constitute an interesting_problem for Physics, determin- 
istic walks in regular 0,131 disordered media 0, H, 0| 
present very interesting results, as an application to for- 
aging [1, [1^ . The memory in random walks has the 
effect of changing the behavior of the gaussian displace- 
ment distribution prj . Here, we are interested in under- 
standing fundamental aspects of a partially self-avoiding 
deterministic walk algorithm, known as the tourist walk 
(TW) [H, [H, [il. These walks, that are described be- 
low, have been applied to characterize thesaurus [l4| . as 
a pattern recognition algorithm [Tsj and image analy- 
sis [Ulil. 

Consider N points (sites, cities) randomly distributed 
inside a c?-dimensional hypercube with unitary edges. 
The distance Dij between any two points Si and Sj is 
calculated via euclidean metrics. The walker leaves a 
given point and moves obeying the deterministic rule of 
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going to the nearest point (shortest euchdean distance), 
which has not been visited in the /i preceding steps. This 
rule produces trajectories with an initial transient part 
of t steps and a cycle of p steps as a final periodic part. 
Once trapped in a cycle, the walker does not visit new 
points anylonger. Short transient times and short period 
cycles limit exploration of the medium by the walker. 

Analytical results [3| could be obtained for (i) mem- 
oryless walkers in the deterministic [l^ and stochas- 
tic 0, [2l[ versions of the TW and for (ii) determinis- 
tic walk with arbitrary memory in one-dimensional sys- 
tems [2^ 1. Here we consider the memory effect in deter- 
ministic walks in a mean field approximation. 

The deterministic TW, with memory /i = 0, is triv- 
ial since the walker does not move at each time step, so 
that the transient-time/period joint distribution is sim- 
ply: SQ'^J{t,p) = St.oSp,!, where 5ij is the Kronecker 
delta. With memory fi = 1, the walker must leave the 
current site at each time step. The joint distribution 
S'f2{t,p) is obtained considering the trajectories of a 
tourist leaving from all sites of a given map and statis- 
tics is performed for different realizations (maps). For 
iV 3> 1, the transient-time/period joint distribution is 
obtained analytically for arbitrary dimensionality [l^: 
S['^^\t,p) = [(t+/-i)r(l+/-i)/r(t+p+/-i)]<Sp,2, where 
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T{z) is the gamma function and Id = /i/4[l/2, {d+ l)/2] 
is the normaUzed incomplete beta function. This case 
does not lead to exploration of the random medium since 
after a short transient, the tourist gets trapped in pairs 
of cities that are mutually nearest neighbors. 

Interesting phenomena occur when the memory values 
are greater or equal to two (/i > 2). In this case, the 
cycle distribution is no longer peaked at Pmin = M + 
1, but presents a whole spectrum of cycles with period 
P > Pmin, with possible power-law decay [13, [3] i which 
favors exploration of the medium by the walker. The 
elucidation of this intringuing broadening of the cycle 
period distribution is our main objective. 

As the medium dimensionality d incresases, the cor- 
relations between the distances Di j become weaker 
and weaker, so that, in the high dimensionality limit 
(d oo), the distances can be considered independent 
random variables, uniformly distributed in the interval 
[0, 1] [13, m, m, [H, [13. This is the mean field model 
named Random Link (RL), where two euclidcan con- 
straints still remain: (i) the distance from a point to itself 
is null, Di^i = 0, and (ii) the forward and backward dis- 
tances are equal, Di j ~ Dj.i- Breaking these constraints 
leads to the Random Map model (RM) [l^. which is a 
mean field approximation for the Kauffman's model [29t . 
The neighborhood statistics for these models have being 
analytically studied in Ref. [sO- 

In this paper, we obtain analytical results for the TW, 
with memory /i = 2 in the d ^ oo medium, i.e. the RL 
approximation. These results enable us to explain the 
main mechanism which makes the fj, = 1 and /i > 2 sit- 
uations so distinct. Also, they permit us to estabilish a 
relationship between the mean fields RL and RM models. 
The walks with memory fi = 2 in the symmetric indepen- 
dent random distance case (RL model) is equivalent to 
memoryless (/x = 0) walks in the assymmetric indepen- 
dent random distance case (RM model), which has been 
already solved in Ref. [l^. Throughout this relationship 
between RL and RM models, we show that the decay for 
the cycle period distribution in the RL model is a power 
law (X p~^. 

Also we are able to explain the reason of the already 
observed numerically abrupt change in the in the tran- 
sient/period joint distribution for null transient t = 0. 

The presentation of these results are briefly skechted 
in the following. In Sec. |TT1 we calculate the probability 
S2^j{n) for the tourist, with memory /i = 2, to visit n 
distinct sites before the first passage to any already vis- 
ited site, walking on the RL model with N sites. We 
start calculating the complementary cumulative distri- 
bution F2^i (h) (upper-tail distribution). Next, through- 
out an analogy to the geometric distribution, we obtain 

the revisit p^^Ji{j) (first passage) and exploration q^^i{j) 
probabilities. Using an alternative derivation, we obtain 
simpler expressions for these probabilities, which lead to 
a closed analytical expression for Plf^J^ (n). In Scc. lIIIl we 
show that the probability for the tourist to be trapped 



into a cycle when revisiting a site is 2/3, which is counter- 
intuitive. This result (combined to previous ones) allows 
us to obtain the complementary cumulative distribution 
^2^1 total number n of visited sites (until the 

walker enters an attractor). In Sec. IIV[ we obtain the 
joint distribution S'^Ji{t,p) of transient time t and cy- 
cle period p and show the drastic difference between the 
t = and t 7^ cases. Final remarks are presented and 
future studies are proposed in Sec. |Vl 



II. DISTRIBUTION FOR THE NUMBER OF 
EXPLORED SITES BEFORE THE FIRST 
PASSAGE (REVISIT) 

Consider that the tourist, who performs a walk with 
memory — 2 on the RL model with N points, has 
visited n>3 = /i+l = distinct sites and then 

revisits one of these sites. Aiming to obtain the distribu- 
tion 82^1 (n) of the number h of sites visited before the 
first passage, we start calculating the complementary cu- 
mulative (upper-tail) distribution 

N 
k—n 

i.e., the probability for the tourist to explore at least n 
distinct sites, before the first revisit. 

In the schema of Fig. [1] the tourist leaves from a given 
site si (first step, j ~ 1) and follows the trajectory si, 
S2, . . . ; Sfj, exploring n = 9 distinct sites, with no revisit. 
For 1 < i < n — 1 , let us denote 

• Xi the distance between the consecutive sites Si 
and Si+i in the trajectory (thick continuous lines 
of Fig. [ID, 

• yi,k the distances between the site Si in the tra- 
jectory and other sites outside the trajectory (thin 
continuous lines of Fig. [1]). 

• Zi^k the distance between the non-consecutive sites 
Si and Sfc in the trajectory (slashed hues of Fig. [1]), 

By definition of the RL model, all these distances Xi, yi^k 
and Zi^k has uniform deviate in the interval [0, 1]. 

The conditions for the tourist to follow the trajectory 
si, S2, . . . , Sft in the first n steps are 

1 . in the case /i = 1 (already solved in Ref. [l9l ) , the 
distances Xi must obey the relation Xn^i < Xn^2 < 
• • • < xi, once the tourist stops exploring new sites 
when Xi^i > Xi, giving rise to a cycle of period 
p = 2. But for the case n = 2 addressed here, each 
distance Xi may vary unrestrictly in the interval 
[0, 1], because the memory fi = 2 forbids the tourist 
to move backward from Si+i to Si (even if Xi^i > 
x^). 
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Figura 1: Schematic representation of a walk with at least 
n = 9 sites visited before the first passage. The walker leaves 
from the site si and follows the trajectory si, S2, S3, . . . , sg. 



2. when the tourist is about to walk the distance Xi 
(and move from Si to Si+i) there exist N — i non- 
explored sites at his/her disposal. 

3. for each site s;. all iV— (?? — 1) distances yi_k must be 
greater than Xi. The probability for this to occur 



IS 



N-n+l 

= (1 - The only ex- 

ception is the site Sfj_i, which has N — n distances 
Un-i^k connected to it (see Fig. [TJ where sa-i cor- 
responds to Sg). 

4. to avoid shortcuts and revisits, each distance Zi^k 
must be greater than both a;,; and Xk- 

These conditions lead to the following chained inte- 
grals: 



n-2 „i 



TT / dx,(7v-j)(i~x,)^-"+' 



dxf,-i{N - h + 1)(1 - Xn-i) 



n n / ^ d..,.. 

i=r fc^i+2-'max(^,,^fc) 



TV- 



(1) 



It is worthwhile to mention that we have made no ap- 
proximation yet, hence Eq. [1] yields exact results even 
for small values of iV, as Tab. |T] shows. 

However, the function max(a;i,a;fc) in the lower limits 
of the integrals in Zi^k makes it difficult to solve Eq. [1] 
once we should consider all possible (fi — 1)! orderings 
of distances Xi. In the following, we will consider the 
thermodynamic limit (N ^1) and make some approxi- 
mations to solve Eq. [TJ 

For a better visualization, notice that the integrals in 
Zi^k refer to the slashed lines of Fig. [TJ Observe that 



27 

32 
9 459 
17 248 
107 301 
R09 fiOO 



0,15625 
0,29534 
0,33785 
0,21056 



standard 
mean error 



0,15624 1 • 10'^ 

0,29535 1 ■ 10"^ 

0,33784 1 ■ 10"^ 

0,21056 1 ■ 10"^ 



dif. (in 
difference std-error) 



7- 10"^ 
2 ■ 10"^ 
1 • 10"^ 
3- 10-*^ 



0,62 
1,13 
0,82 
0,22 



Tabela I: Numerical validation of Eq.[TJ The columns f'2^ri(^) 

and s!fli{n) refer to analitical values and the columns mean 
and standard-error came from numeric simulation. Walks 
were performed on 300 000 000 maps with A'' = 6 points each. 



exactly n — A slashed lines leave from each site, except 
the sites si and sa-i, where n — 3 slashed lines leave 
from, due to the additional distance zi^s-i (thick slashed 
line in Fig. [TJ). To obtain a more regular expression, 
we can eliminate the integral in zi^n-i in Eq. [TJ (with- 
out any harm) and then each variable Zi_k appears ex- 
actly n — 4 times. To justify this elimination notice that, 
due to the deterministic rule of TW, each distance Xi 
is the minimum of — 2 random variables uniformly 
distributed in the interval [0,1]. Therefore, its pdf is 
given by [l^: g{xi) = (TV — 2)(1 — Xi)'^~^ and its mean 
and standard deviation are: Xi = 1/(TV — 1) w 1/A^ and 
cr^r. = ^/iN ~2)/[N{N -1)^] « 1/iV, so that, in the 
limit N ^ 1, Xi assumes values close to and the value 
of the integral /j^ax(a:i x- 1) d^^i.n-i is typically close to 
1. 

Changing the exponent of xa-i from N — n to N — 
fi -|- 1, all the variables Xi are raised to the same power. 
The resulting expression is algebraically symmetric with 
respect to the variables xi, X2, ■ ■ ■ , Xn-i, what means 
that all possible {n — 1)! orderings occur with the same 
probability. Thus, one can consider the specific ordering 
xi < X2 < ■ ■ ■ < Xn-i and rewrite Eq. [TJ without using 
the inconvenient function max(): 



(n-1)! 



T\{N-l) d.Tl(l-.Tl)^-"+l 

n / dx,{i~x,r-^+'-^ 

/ d.Tft_l(l-a;ft-l)^-^ 



(2) 



where we emphasize that the extra factor {fi — 1)! takes 
into account all possible orderings of the variables Xi. 

The exponent of xi may be changed from A'^ — n -I- 1 
to N — n aiming the exponents of xi, X2, ■ ■ ■ , Xn^2 to be 
in an arithmetic series. One then calculates the integrals 
of Eq. [2Jto have at last: 



{n - 1)1{N - 1){N - 2)(Af - 3) • • • 
{N - 2){2N - 4)(3iV - 7)(4]V - 11) ■ 
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{{n^l)N -[{fi-l)n/2 + l]} 



n 

h 

n 

J=4 



k{N -k) 



kN -k{k + l)/2~l 



N-j + 1 



N-j/2^1/{j-l) 



(3) 



where we have called j = k + 1 and the lower limit of the 
productory was changed from j = 2 to j = 4 because the 
factors for j = 2 and j = 3 are physically meaningless, 
as we shall argue in the Subsec. Ill Al 

Finally, the distribution of n is calculated from the one 
step difference of the upper-tail distribution: 



1 - 



N-fi 



n 

J=4 



N - (n+ l)/2 - 1/n 

N-j + 1 
N-j/2-l/{j-l) ■ 



(4) 



The expression of Eq. [3] is similiar to the one obained 
for /i = 1 (using Eqs. 9 and 10 of Ref. [Tp and calling 

h = t + 2): F(^)(n) = [n;=3W^]/(^-l)!- The main 
difference is the presence of the factor 1 / (n — 1 ) ! , because, 
for /i = 1 , one must to consider only the specific ordering 
< Xn-2 < ■ ■■ < Xi- 
At this point we arc able to understand the major role 
played by the memory in this partially self avoiding walk. 
For /.i = 1, the walker must go to the nearest neighbor. 
The extremal statistics is behind this dynamics. But, for 
instance, forbidding the walker to return to the last vis- 
ited site, this opens up the possibility to go to the first or 
second nearest neighbor, which transforms the extremal 
statistics to the combinatorial statistics. Mathematically, 
this is expressed by the absence of (ri — 1)! in Eq. [31 



factor explains the abrupt change in the exploratory be- 
havior between ^ = 1 and ^ = 2 cases: on one hand, for 
/i = 1 the exploration probability (in the thermodynamic 
limit) decreases harmonically along the trajectory; on the 
other hand, for /i = 2 the exploration probability tends 
to 1 when N —^ oo. 

Once the memory fi = 2 assures the tourist to explore 
at least fimin = fi+1 = 3 sites, it only makes sense to de- 
fine exploration probability from the 4th step. In fact, for 
the first step (j = 1) Eq.Odoes not have a defined value, 
for the second step it yields qi1l{2) = {N -l)/{N-2) > 

1, which is an absurd, and for the third step 0^2^; (3) = 1- 
To take into account the proper physical content, we pre- 
viously changed lower limit of the products of Eq.[3]from 
j ~ 2 to j ~ 4. Its interesting to mention that for the 
step j — N + 1 (after the tourist explores all N sites), 
Eq. [5] correctly yields qi^JiN + 1) = 0. 

Since in the jth step there are j — 3 sites equally prob- 
able to be revisited and P2^](j) is the probabihty for 
the tourist to revisit any one of these sites, in the limit 
N ^ j 1 the probability Pri for the tourist to revisit 
a specific site Sk is 



1 



Prl 



3-3 
1 

2N ' 



1 



j/2 - 1 - l/jj - 1) 
j-3 N-j/2-l/(j-l) 



(6) 



which is half the probability for he/she to explore a spe- 
cific new site [namely (jri — 1/{N — j) w 



B. Alternative Derivation 

In the following wc obtain for iV ^ 1 simpler expres- 
sions for the first passage and exploration probabilities, 
via an alternative reasoning. From these probabilities, 
we obtain closed analytical expressions for F!fii{n). 



A. Analogy to the geometric distribution 

Making an analogy to the geometric distribution, we 

can write Eq. g] as S^^J, (n) = p%\ {n + 1) UU C ) 
where 



N-j + l 



N -j/2-l/{3-l) 



the 



(5) 



and 



Co-) 



exploration probability in the jth step 

= 1 — 92^; (j) is the revisit probability in the 
jth step. We remark that the expression of Eq.[5]is simi- 
lar to the one obtained for /i = 1 (adapting Eqs. 9 and 10 
of Ref. [l^ from their original concept of subsistence prob- 
ability to the concept of exploration probability handled 
here): {j - l)q[l\{j) ^ {N - j + l)l{N ^ ] /2). The main 
difference is the extra factor j — 1 , which is a consequence 
of the restriction Xn-i < < • • ■ < xi. This extra 



1. First Passage and Exploration Probabilities 

Supose that the tourist has traveled along the trajec- 
tory si, S2, Sn {n > 3) without any revisit. Let 
us first reobtain the probability Pri for the tourist to re- 
visit a specific site Sfc (outside the exclusion window, i.e., 
fc < 71 — 2) in the following step. To do this, consider the 
following constraints (see Fig. [1]): 

1. the distance Zn^k must be smaller than x^. 

2. once in the (fc + l)th step, the tourist came from 
site Sk to Sfe+i, the distance za^k is greater than the 
distance Xk- 



In brief, z^^k 

Xk ^ '^h.k ^ <C 1 



must vary between Xk and Xn^ so that, < 
C^i < 1. 



5 



Once the pdf of each distance Xi is g{xi) = {N — 
2)(1 — Xi)'^~^ and Zn^k has uniform deviate (by defini- 
tion of RL model), for iV 3> 1 the probability Pri is 
given by: pri = P{xk < z^^k < Xn) ^ Jq dxk{N - 
2)(1 - ^^n{N - 2)(1 - dzr,,k = 



1){2N 



3)] « 1/(2A^), which agrees to 



{N ~ 2)/[{N 

Eq.m 

For a generic step j there are j — 3 sites susceptible 
to be revisited so that the first passage and exploration 
probabilities for this step are: pi^Jiij) = (j — 3)/(2iV) = 

1 — 92^;(j), which is an approximation for Eq. [5l leading 
to 



r(27V) 

r(2iV-n + 3)(2Ar)"-3 
which is a closed analytical form for Eq. [3l 



2N 



(7) 



2. Exponential Form (Cumulative Half Gaussian) 

In the limit TV 3> 1, the exploration probability may 
be written as qf^\{j) = [1 - 1/{2N)]^-'^, so that Eq. [7] 
assumes its exponential form 



ncio)^(i-5^)' 



where 



i=4 



3) 



-[(«-3)V(4W)][l + l/(fi-3)] 



(n - 2)(n- 3) 



,(8) 



(9) 



has a simple physical interpretation. It is just the num- 
ber of distances Zi^k between non-consecutive sites of tra- 
jectory. Notice that the trajectory of Fig. [T] is topolog- 
ically equivalent to a (n — l)-sided polygon, which has 
(ri — l)(n — 4)/2 diagonals. All these diagonals plus the 
side siSfi-i totalize Q = {n — 2){h — 3)/2 paths (slashed 
lines of Fig. [T]), which allow revisit. 

For n— 3 ^ 1, one can disregard 1/(?t-— 3) in Eq.[51 lead- 
ing to a half gaussian: y = i^2^i^(n) = e"!'""^^/^'^!^/^, 

indicating that the scaled variable is x = (n—nmin) / 
with fimin = /-t + 1 = 3, leading to the universal curve 
y = with x >0. We only have kept fimin to com- 

pare to a possible generalization of these calculations for 
the case of short memory /i <C A^. 



III. DISTRIBUTION OF THE TOTAL NUMBER 
OF EXPLORED SITES 

Up to this point we were focused on the number h of 
sites explored before the first revisit. In the TW with 



11 = 1, the revisit implies the tourist has entered an at- 
tractor of period p = 2 [l^, but with n = 2, the revisit 
does not implies capture. In what follows we calculate 
the probability pt for the tourist to be trapped during 
a revisit and then obtain the capture p'^Jiij) and sub- 
sistence 92^; (j) probabilities and the upper-tail distri- 
bution F2^^ (n) for the number n of sites visited in the 
whole walk. 



A. Trapping Probability 

Let us recall Fig. [1] and consider that the tourist has 
traveled along the trajectory si, S2, • • • , Sfi without any 
revisit. Assume that in the following step he/she revisits 
site Sk (outside the memory window, k < h — 2). Due 
to the deterministic rule, two situations may occur: (i) 
if Xk < Xk-i, the tourist moves forward to site Sk+i and 
is trapped by an attractor of period p = n — fc + 1; (ii) 
if Xk-i < Xk, the tourist moves backward to site Sk-i 
and escapes from the attractor. Therefore, the walker 
trapping or escaping depends on which distance Xk~i or 
Xk is shorter. The only exception is a revisit to si, when 
the tourist is unconditionally trapped, leading to a tra- 
jectory with a null transient time {t = 0) and a cycle of 
period p = h. 

Taking into consideration that all (n — 1)! possible or- 
derings of the distances xi, X2, xn-i are equally 
probable, one could naively conclude that the trap- 
ping probability would be pt = P{xk < Xk-i) = 1/2. 
Nonetheless, numerical simulations of this system have 
refuted this expectation, pointing out that this probabil- 
ity is in fact pt = 2/3. 

To understand this result, we first show that the prob- 
ability Pvif) for the tourist to revisit a specific site Sk is 
proportional to the rank r occuped by the associated dis- 
tance Xk (between sites Sk and Sk+i) when one reorders 
the distances xi, X2, . . . , Xn-2 decrcasingly (so that Xk is 
the ?'th greatest one) . Secondly, we show that the proba- 
bility Pt{r) for the tourist to be trapped when revisiting 
the site Sk is proportional to r — 1. Finally, from Pvir) 
and Pt{r) we prove that pt = 2/3. 



1. Order Statistics 

Let us recall some tool about Order Statistics. Given 
a sample of M variates Xi, X2, . . . , Xm, reorder them 
so that A(i) > A(2) > . . . > A(jv/) . If X has pdf 
g{x) and cumulative distribution G{x) = J^^dx'g{x'), 
then the pdf Kix) of A(,,) is hr{x) = M\[G{x)]^^-''[l - 
G{x)Y-^gix)/[ir - 1)!(A/ - r)!], for r = 1, 2, . . . , M. 

Resuming the TW with /j, = 2 on the RL model, 
each distance Xi has pdf given by g(x) = (A — 2)(1 — 
then its cumulative distribution is Gix) = 

dx' g(x') = 1 — (1 — x)^~^ and the pdf of xj^) 
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is hr{x) = n! [l-{l-x)^-^y 
2)(1 -a;)^-3/[(^ -!)!("- 0!]- 



[(1 - x)^-^Y-\N- 



2. Rank-revisit and Rank-trapping Probabilities 

Again, consider that the tourist has traveled along 
the trajectory si, S2, ■ • ■ , sa (without any revisit). Let 
us calculate the probability Pv{r) for he/she to revisit 
the site = (with associated distance X(^r) = Xk) in 
the next step. Once is the nearest site, the distance 
Zfi has pdf given by g{x) = {N — 2)(1 — x)^^^ and 
once the tourist came from site S(r) = Sk to Sk+i in the 
(fc + l)th step, the distance Zn^(^r) is certainly greater than 
Thus Py{r) (X P(2;ft.(r) > = h\/[{r - l)!(n - 

r)!] dx [1 - (1 - x)^-'] [(1 - xr-^] '-\n 
2)(1 - x)^-^ dz(iV- 2)(1 - z)^-^. Evaluating the in- 
tegral in z and calling y = (1 — x)^^^ the above equation 
is rewritten as: Pv{f) oc n!/[(r— l)!(ri — r)!]B(?T, — r + 1, ?■+ 
1) = n!/[(r - l)!(n - r)!](n - r)\r\/{n + 1)! = r/(n + 1). 
This expression is not the probability P„ itself. Instead, 
it only gives the dependence of on r. 

Normalizing P^ over 1 < r < ri — 2, one has 

where n — 2 is the number of sites available to revisit 
(the sites Sa and Sa-i are forbidden by memory) and the 
normalization factor "Y^^Ji j = {n — 2){n— l)/2 is simply 
the sum of all n — 2 ranks. 

The result of Eq. [TO] does not contradicts Eq. [51 since 
Eq.[6| gives an approximated probability for the tourist to 
revisit a specific site Sk , regardless its associated distance 
Xk = a^(r)j while Eg. 1101 gives the conditional probability 
for the tourist to "choose" the r-ranked site S(r) during 
a revisit after exploring fi distinct sites. 

Once the tourist had revisited site Sk (or equivalently 
S(r)), the probability Pt(r) for he/she to be trapped also 
depends on the rank r. The trapping condition is that 



Xk-i must be greater than Xk- Since Xk 



is the rth 



greater distance, there are only r — 1 remaining distances 
(among ri — 3 ones) greater than Xk- Thus, 

■ - 1 



Pt{r) 



n — 3 



(11) 



Combining Eqs. [TO] and [TTl the probability for the 
tourist to be trapped when visiting a specific site S(r) 
is: P^,(r)Pt(r) = 2r(r-l)/[(n-l)(n-2)(n-3)]. Thus, 
the probability for the tourist to be trapped when revisit- 
ing any site is pt = Yl^Zi Pv{''')Pt{r). Calling m = n — 2 
and evaluating ''(^ ~ 1) = m{m? — l)/3 one finds 

the trapping probability 



Pt = 2/3 . 



(12) 



We remark that this result has been obtained without 
any approximation, and numerical simulations agree to 
it even for small values of N . 



B. Capture and Subsistence Probabilities 

Combining the probability pri for the tourist to revisit 
a specific site Sk (Eq. [S|) and the trapping probability pt 
(Eq. [TO)) , one obtains the probability pri for the tourist 
to revisit Sk and be trapped: 



Pri = PrlPt 



1 2 



1 



2A^3 3iV 



(13) 



Since in the jth step there are j — 3 sites available to 
revisit, the capture (i.e., revisiting any site and being 
trapped) and subsistence (i.e., exploring any new site or 
revisiting any site and not being trapped) probabilities 

in the jth step are: p'i'^\{j) = [j - 3)/(3^) = 1 - gg](j) 
and the upper-tail distribution for the number n of sites 
explored by the tourist in the whole trajectory is 



r(3iv) 



(14) 



r(37V- n-h3)(3iV)«-3 ' 
which is analogous to Eq. [7| 

1. Comparison to RM model with /j. = 

The expression of Eq. [T4| is similar to the one ob- 
tained for the RM model with memory ^ = p^ : 

-Fo,™(") = T{N)/[r{N - n)N"]. This result explains 
the non-trivial equivalence observed between RL model 
with N points and memory fi = 2 (memory effect) and 
RM model with 3iV points and memory fi — (effect of 
distance symmetry break), when one compares the dis- 
tributions for the total number n of sites explored by the 
tourist. 

Notice that, taking both models with points each, 
in RL with /i = 2, at each step, the probability for 
the turist to revisit a specific site and be trapped is 
Pri w 1/(3A^); and in RM with ^ = Q, this probabil- 
ity is Prm = l/N. Therefore, taking RL with N points 
and RM with 3A^ points equals these probabilities and 
justifies the equivalence. 

2. Exponential Form 

In the limit N ^ 1, the subsistence probability is 
rewritten as 92^;(j) = [1 — 1/{3N)Y~^ and one ob- 
tains the exponential form of Eq. [TH namely F2^i (n) = 



nu&) = [1-1/(3^)]^ 



(n- 2)(n- 3)/2. 
Rather than differentiating distribution 

S^^j (n) for the number n of sites explored in the whole 
trajectory is more precisely obtained by imposing the 
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tourist to explore n distinct sites and then be captured 
in the next step (i.e., revisit any site and be trapped): 



n — 2 t7i-2)(Ti-3)/2 
— g 3JV 

3iV 



(15) 



For n ^ I, calhng y = \J SNS^^j (n) and x = (n — 

n„iin) /VSiV (with rimin — IJ, + I = 3) one obtains the 
universal plot for this system: 



y 



(16) 



with cc > and mth moment (cc™) = 2'"/2r(m/2 + 1), 
where we see that normalization is assured by (x*^) — 
1. The mean value is (x) ~ A/7r/2 and the variance 
(x^) - = 2 - 7r/2. Fig. [2] exhibits a plot of Eq. [H 
and experimental data. From this figure, or calculating 
analytically, one obtains that the mode is unitary. 



▲ 


N = 10^ 


■ 


N = 10^ 


•*- 


TV = 10^ 


• 


TV = 10" 




- Theoretical 




Figura 2: Finite size effect for the distributions and con- 
vergence to the universal curve: y = x e ^ ' , with x — 

€1< 



(n - nrmn)/V3N and y = VMS^^Hin). 



IV. TRANSIENT AND PERIOD JOINT 
DISTRIBUTION 



(t + p-2)(t + p-3)/2 

e ™ (17) 



site Si]: 

where Sij is the Kroncckcr delta. Fig. [3] exhibits a plot 
of Eq. [l7]for TV = 1000 points. 



^1^1 (^,P) 




80 100 ^0 



Figura 3: Transient and period joint distribution for a map 
with TV = 1000 points in the RL model with fi = 2. 



A. Transient time marginal distribution 



The transient time distribution is calculated sum- 



ming Eq. [17] over all possible periods, i.e. ^2^1(1) 



J2p=3 ^2 rii^^P)- limit TV ^ 1, this summation 

can be approximated by the integral 



5/2 
1 



dpC;(i,p) 



— I /— crfc 
2 J \ 6n'^^^ 



/6TV 



where the lower limit 5/2 is due to a Yates continuity cor- 
rection (which other than improve the integral approxi- 
mation, make the analytical form quite simpler) and the 
upper limit has been extended to infinity to make calcula- 
tion easier (with no harm, because Eq. [TT] yields despica- 
ble values for p > TV) and erfc(x) = (2/^/7r) dxe~^ 
is the complementary error function. 



The transient /period joint distribution 5*2^; (t,p) can 
be obtained similarly to Eq. 1151 by imposing the tourist 
to expore n distinct sites and then rcvist the specific 
site Sfe (instead of any site) and be trapped, giving rise 
to a tracjetory with transient t = k — 1 and period p — 
n — k+1. We notice that the relevant variable is t+p = n. 
Hence, s'^Ji{t,p) is obtained multiplying F2^j[t~\-p) by 
Pri (Eq. [131) [or by pri (Eq. ^ in the case t = 0, since 
the tourist is unconditionally captured when rcvisting the 



B. Cycle period marginal distribution 

Similarly, the period distribution is 

JV-3 



s'2>) = E€ht^p)= / '^is^^'tlit.p) 

t=o -^-1 
p-7/2 




pV(6JV) 
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where the lower limit —1 is due to both Yates continuity 
correction and a compensation for the half extra degree 
in t = 0. The mean period value is p = y^SirN/S and 
standard deviation is dp = ^^(2 — 37r/8)A^. For p <C 
VGN, the decay follows a power law S2^]{p) (xp~^. 

V. CONCLUSION 

In this paper, we have analytically obtained the statis- 
tical distributions for the deterministic tourist walk with 
memory fi = 2 on the random link model. 

The distribution for the number of sites explored before 
the first passage has been compared to the one previously 
obtained for the case /x = 1, elucidating the mechanism 
that strongly increases the tourist's exploratory behav- 
ior. On one hand, for = 1 the distances travelled at 
each step must obey xi > X2 > . . ., leading to a local- 
ized exploration. In the thermodynamic limit, the mean 
number of explored sites is then n = e = 2.71828 . . . and 
the exploration probability decreases harmonically along 
the trajectory. This dynamics is due to the underlying 
extremal statistics. On the other hand, for fi = 2 the 
distances xi,X2, ■ ■ ■ are unconstrained, leading to an ex- 
tended exploration: n is proportional to N^^"^ and the 
exploration probability tends to 1, when N oo. This 
dynamics is due to the underlying combinatorial statis- 
tics. The factor {n — 1)1 in Eq. [2] represents the change 
from the extremal statistics to the combinatorial one, 
which makes the Sp 2 distribution of ^ = 1 to broaden to 
a wide (l/p) distribution for /i > 2. 

Throught the trapping probability pt = 2/3 (which 
value is counterintuitive), we obtained the capture and 
subsistence probabilities and a closed form to the comple- 
mentary cumulative distribution for the number of sites 
explored in the whole trajectory. This distribution is 
analogous to the one obtained for the random map model 
with /.t — 0. This result explains the equivalente between 



these mean field models (RL with N points and memory 
/i = 2; and RM with 3A^ points and memory /i = 0). 

For a large number of sites {N ^ 1) in the random 
medium, the distribution S2^l{n) of having n distinct 
sites visited by the tourist with memory ^ = 2 in the 
random link model is universal y = xe~^ ' with y = 
^/3iV5^^;V) and x = {n - 3)/VM. 

The transient time t and cycle period p joint distribu- 

tion Si'^]{t,p) = e[(*+P-3)'/(3A^)]/2/[7v(3-5t,o)] has been 
obtained noticing that the revclant variable is approx- 
imatively given hy t + p = n. The marginal distribu- 
tions are also universal. For the transient time one has: 
y = [1 + (5(x)/2]erfc(x) with y = y^m/^s!^^J{t) and 

X = t/ \J 6N and for the period distribution: y = erfc(x), 
with y = y/m/^S^^^J{p) and x = {p - 7/2)/VQN. We 
have shown that the discrepance in the null transient time 
distribution {t ~ 0), when compared to the subsequent 
ones {t > 0), is due to the higher capture probability the 
starting site si has [namely, Pri = 1/(2A'^)] when com- 
pared to the others else [pri = 1/(3A^)]. Wc also have 
shown that the period distribution decays according a 
power law 52^; (p) oc p~^. 

Future studies concern the consideration of higher 
memory values in the random link model and the under- 
standing of the connection with the random map model. 
As the memory increases, we expect a transition from the 
closed periods to non-closed ones (chaotic phase). We are 
interested in understanding the role of finite dimension- 
ality of the system. 
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